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We report on a new constraint on gravitylike short-range forces, in which the interaction charge 
is mass, obtained by measuring the angular distribution of 5 A neutrons scattering off atomic xenon 
gas. Around 10 7 scattering events were collected at the 40 m small angle neutron scattering beam 
line located at the HANARO research reactor of the Korean Atomic Energy Research Institute. The 
extracted coupling strengths of new forces in the Yukawa-type parametrization are g 2 = (0.2 ± 6.8 ± 
2.0)xl0 -15 GeV -2 and g 2 = (—5.3±9.0 jf 2 g)xlO -17 GeV -2 for interaction ranges of 0.1 and 1.0 nm, 
respectively. These strengths correspond to 95% confidence level limits of g 2 < (1.4 ± 0.2) x 10 -14 
GeV -2 and g 2 < (1.3 ± 0.2) x 10” 16 GeV -2 , improving the current limits for interaction ranges 
between 4 and 0.04 nm by a factor of up to 10. 


Extensions of the standard model of particle physics 
have long been discussed. Some of these theories, based 
on supersymmetry or extra space dimensions, naturally 
involve gravity or gravitylike interactions even at low en¬ 
ergies. Several models predict new bosons mediating 
gravitylike forces that couple to mass, baryon number, 
or, in the case of grand unification models, to the differ¬ 
ence of baryon and lepton numbers [H-@|. This class of 
models induces modifications to the Newtonian inverse- 
square law of gravitational interactions, and may also 
cause violation of the weak equivalence principle. These 
proposals motivate us to search for new gravitylike forces. 
Comprehensive reviews of such theories, and the forces 
they predict, can be found in Refs. Bit 

The forces due to such new bosons can be simply mod¬ 
eled by a Yukawa-type scattering potential written in 
natural units as 
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where g 2 is a coupling strength, Qi are coupling charges, 
and /r is the mass of the boson mediating the force. Such 
models can be considered in the g 2 — /r or g 2 — A param¬ 
eter space, where A s 1 /pi is the interaction range. Cur¬ 
rent experimental limits at 95% confidence level (C.L.) 
for interactions that couple to mass are shown in Fig.|T| 
Constraints A and B were obtained by microscopic ex¬ 
periments that precisely measured interactions between 
neutrons and atoms 10j,[ll|. The experimental method 


used to achieve the results reported in this Letter follows 
a similar approach. Constraints C to I were obtained 
by macroscopic tests, searching for non-Newtonian forces 
between test masses using techniques such as torsion bal¬ 


ances and microcantilevers 12}-1 


This experiment was performed at the 40 m small angle 
neutron scattering (SANS) beam line located at the HA¬ 
NARO research reactor of the Korean Atomic Energy Re¬ 
search Institute [Til }. Figure[2] shows a schematic drawing 
of the experimental apparatus. The angular distribution 
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FIG. 1. Experimental limits at 95% C.L. for new gravitylike 
forces in the Yukawa-type parametrization space. Constraints 
A and B were obtained by measuring neutron-atom scattering 
0,0, while constraints C to I were obtained using macro¬ 
scopic methods 0-Q3|. Theoretical expectations from an ex¬ 
tra (7(1) gauge boson at different symmetry breaking scales 
Aj/(i) (~ 246 GeV and 10 TeV) Q,i , and from a baryon num¬ 
ber gauge field in the bulk of extra space dimensions aa are 
shown as dashed lines and the hatched area, respectively. 


of neutrons scattered by a xenon gas target was mea¬ 
sured, and deviations from the expectations from known 
interactions were used to set limits on additional, un¬ 
known, interactions. 

Neutrons with an average wavelength of 5 A and a 
spread of 12% (FWHM) were used, selected by appro¬ 
priately setting the rotation rate and tilt angle of a heli¬ 
cal slot velocity selector installed at the entrance to the 
beam line. The neutron beam intensity was monitored 
by a 3 He-£lled sampling gas counter located immediately 
downstream of the velocity selector. Two collimators 
made of sintered B4C with a circular aperture of diame¬ 
ter 22 mm were placed 12.2 m and 20.0 m downstream of 
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FIG. 2. A schematic view of the experimental setup (not to 
scale). Neutrons of wavelength 5 A and a 12% spread were 
chosen by a velocity selector. Two collimators of diameter 
22 mm produced a neutron beam with 3 mrad. divergence, 
which impinged on a xenon-filled target chamber of length 
250 mm. Scattered neutrons were measured by a MWPC 
with pixel size 5.1 x 5.1 mm 2 located 3.11 m downstream of 
the gas target. The xenon gas was purified to better than 10 
ppm by a getter system using a carefully controlled flow rate. 


the gas counter, defining the beam divergence of 3 mrad. 


A cylindrical target chamber with an inner volume 
of length 250 mm and diameter 130 mm was placed in 
the sample chamber immediately downstream of the sec¬ 
ond collimator. The target chamber was connected to a 
getter-based purifier which reduced contaminants in the 
xenon gas to less than 10 ppm, a vacuum system consist¬ 
ing of turbomolecular and scroll pumps, and xenon gas 
bottles. The neutron entrance and exit windows of the 
target chamber were made of a single Si crystal [20j, and 
only metals were used in the vacuum seals to avoid con¬ 
tamination by outgassing. The chamber design allowed 
the use of a maximum xenon gas pressure of 2 atm. 

Neutrons scattered by the xenon gas target were mea¬ 
sured by a multiwire proportional chamber (MWPC) 
containing a mixture of 60% 3 He and 40% CF 4 . The 
active volume of the MWPC was 980 mm (horizontal) 
x 980 mm (vertical) x 63.5 mm (depth), and its de¬ 
tection efficiency was 80% for neutrons of wavelength 5 
A. The MWPC consisted of a single anode wire plane 
sandwiched by two cathode wire planes at a distance of 
6.4 mm. Each plane consisted of 192 wires at a pitch of 
5.1 mm, giving a spatial resolution of 5 mm (FWHM). 
The front (rear) cathode plane was oriented vertically 
(horizontally), so the detector can be treated as a two- 
dimensional pixelated detector with pixel size 5.1 x 5.1 
mm 2 . Calibration of the relative detection efficiencies of 
each pixel to an accuracy of 0.5% was performed using 
a standard water sample 21]. The detector was placed 
3.11 m downstream of the center of the target chamber. 
The entire neutron beam path was evacuated to 0.1 Pa 


to minimize scattering off residual gas. 

The differential scattering cross section of unpolarized 
neutrons on xenon atoms due to known and new interac¬ 
tions can be expressed as [ 22 ] 
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where g is the momentum transfer, b c (~ 5 fm) is the 
coherent scattering length, b s (q) (~ 10 -3 fm) is the 
Schwinger scattering length, bi is the incoherent scat¬ 
tering length (negligibly small in the case of the xenon 
atom), and bp (~ 10 -3 fm) is the Foldy scattering length. 
Xem represents the strength of the known non-constant 
term in the coherent scattering processes due to elec¬ 
tromagnetic interactions between the neutron’s spin or 
charge distribution and atomic fields. It is defined as 



FIG. 3. (a) Simulated corrected differential cross sections 

for the constant term h c (9), the electromagnetic interaction 
term h errl (9), and the new interaction term h n ew(9; A = lnm). 
They are normalized to J s h(9)d9 = 1 over the signal accep¬ 
tance S of 16.1 mrad < 9 < 160.8 mrad. Within this sig¬ 
nal region, the distribution due to the new interaction term 
for X ~ 1 nm is clearly different from the known interac¬ 
tion terms, (b) Measured distributions for the xenon-filled 
\gna.m{0)] and empty \g e m P (9)\ chamber. The contribution 
from beam background gb g {9), not shown in the figure, has 
a flat distribution, with around 3 counts/bin. (c) Residu¬ 
als from the reference distribution due to known interactions. 
The data, shown as crosses, are consistent with the case of no 
new forces (y 2 /ndf = 1.4). 
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Xem = {bp + bi)Z/b c , where bj (~ 10 “ 3 fm) is the intrin¬ 
sic neutron-electron scattering length and Z is the atomic 
number. Xem is of order 10“ 2 . The atomic form factor 
can be described to an accuracy of 10“ 4 by the empir¬ 
ical formula /(g) = [1 + 3(g/go) 2 ] -0 ' 5 , where go = 6.86 
A -1 22 1. The additional scattering length due to a new 
Yukawa-type scattering potential is given in the third 
term, where Xnew = m n g 2 QiQ 2 / 2 nb c /j, 2 . With Q i, Q 2 
equal to the neutron and xenon nuclear masses, /i = 200 
eV and g 2 = 10 “ 16 GeV “ 2 give Xnew = 2 x 10“ 3 . Equa¬ 
tion ([3]) is obtained by neglecting terms smaller than 10 “ 4 
fm 2 . 

The expected angular scattering distribution to be 
measured was derived from this differential cross section 
convoluted with the finite beam size, the length of the 
scattering chamber, and the thermal motion of the xenon 
gas, which follows the Maxwell-Boltzmann distribution 
at 293 K. These effects were simulated using the Monte- 
Carlo method, and a corrected differential cross section 
dd/dfl(9) corresponding to this experimental setup was 
expressed as the sum of three terms h c (9), h em (9), and 
hnew(9; /i), corresponding to the constant, electromag¬ 
netic, and new interaction terms of Eq. ([5]) : 


—(9) = N{(l-a*)(l-p)h c (6) 

+a*(l - /3)h em (9) ± /3h new (9; n)} . (4) 

The distributions are normalized to J s dd/d£l(9)d9 = N 
and f s h{9)d9 = 1 over the signal acceptance S of 
16.1 mrad < 9 < 160.8 mrad, and a * : 1 — a* is 
the ratio of scattering contributions due to the electro¬ 
magnetic and constant terms, a* is determined to be 
(1.09 ± 0.01) x 10 “ 4 from measurements of the neutron 
mean-square charge radius (r 2 ) = —0.1161 ± 0.0022 fm 2 
[13. /? is the quantity to be estimated, the fraction due 
to new forces. FigureH a ) shows simulated distributions 
for these three terms. 

The distribution of neutrons was measured for beam 
backgrounds gbg{9) for 3 h, with an empty target cham¬ 
ber g e m P {9) for 36 h, and with a xenon-filled chamber 
gsam(9) for 72 h, separately. The corresponding number 
of counts registered by the flux monitor are Mb g , M emp , 
and M sam . The scattering contribution due to the xenon 
gas sample g{9) was estimated as 

g{0) = g S am{9) - 'l 1 f r Sam gemp{9) 

-M-emp 

( 5 ) 

Mbg 

where 7 is the transmission of the xenon gas, measured 
to be 0.904 ± 0.004. Fig.Hb) shows the distributions 
gsam{9) and g e mp{9). The net distribution from the 
xenon gas is well described by a linear sum of only the 
h c (9) and h em (9) terms (corresponding to /? = 0, or no 
new forces), for which the y 2 per degree of freedom is 



FIG. 4. (a) Evolution of the best fit values g 2 as a function of 
the fraction of collected events k for X = 1 nm. The full data 
set corresponds to k = 1. (b) Evolution of the estimated g 2 
uncertainty <t§ 2 (k). It is well described by the simple model 
CTg 2 (>v) = (f/yTc) © u, with t = 7.6 x 10“ 17 GeV“ 2 and u = 
4.8 x 10 —1 ' GeV“ 2 , shown as the dashed curve. 
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FIG. 5. Comparison of the 95% C.L. limits presented in this 
Letter with those given in QM1. The results improve the 
previous constraints by a factor of up to 10 for interaction 
ranges between 4 and 0.04 nm. The dashed line shows the 
theoretical prediction due to extra U(l) gauge bosons with 
symmetry braking scales of ~ 246 GeV and ~ 1 TeV 


127/90 = 1.4. The residuals from the reference distribu¬ 
tion due to known interactions are shown in Fig. He). 

The measured scattering distribution was fitted by 
Eq. (H for several hypotheses of a new boson’s inter¬ 
action range. As examples, the best fit values of /3 
were /3 = (0.1 ± 2.7) x 10 “ 2 at A = 0.1 nm and /3 = 
(—0.7±1.2) x 10 ” 3 at A = 1.0 nm, corresponding to g 2 = 
(0.2 ± 6 . 8 ) x 10 “ 15 GeV -2 and g 2 = (-5.3 ±9.0) x 10 “ 17 
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GeV^ 2 , respectively. The systematic effect on the ex¬ 
tracted limits due to the uncertainty of the 7 mea¬ 
surement is estimated to be ±15%(A = 0.1 nm) and 
±18%(A = 1 nm), and less than 0.1% due to the uncer¬ 
tainty on a*. Additional effects due to the neutron wave¬ 
length determination, gas temperature, pixel efficiency 
calibration, and atomic form factor modeling were tested 
using a pseudoexperiment technique and confirmed to be 
significantly smaller than the obtained sensitivities. To 
check for unexpected time-dependent systematic effects, 
the evolution with the increasing data set size of g 2 ex¬ 
tracted for A = 1 nm and its error cqp ( k ) are shown in 
Figs.|U)a) andQJb). The errors are well described by the 
simple model ) = ( t/y/n ) © u, where n is the data 
set fraction, and t and u are constants, determined to be 
t = 7.6 x 10 -17 GeV -2 and u = 4.8 x 10 -17 GeV -2 . No 
additional sources of significant systematic uncertainties 
were identified. 

Limits on g 2 at 95% confidence level were evaluated 
using the Feldman and Cousins approach [24]. The ob¬ 
tained upper limit curve is shown in Fig. [5] These results 
improve previous constraints for gravitylike forces in the 
4 to 0.04 nm range by a factor of up to 10. 
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